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>- ■ Recent observations of distant supernovae imply, in defiance of expectations, that the uni- 

^ ■ verse growth is accelerating, contrary to what has always been assumed that the expansion 

^ ■ is slowing down due to gravity. In this paper a general-relativistic cosmological theory that 

o ■ 

O ■ gives a direct relationship between distances and redshifts in an expanding universe is pre- 

Q^' sented. The theory is actually a generalization of Hubble's law taking gravity into account 
O 

d ■ have three phases of expansion, decelerating, constant and accelerating, but it is shown that 

•'-^ ■ at present the first two cases are excluded, although in the past it had experienced them. 
>< 

■ Our theory shows that the universe now is definitely in a stage of accelerating expansion, 

confirming the recent experimental results. 



by means of Einstein's theory of general relativity. The theory predicts that the universe can 
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The general-relativistic theory of cosmology started in 1922 with the remarkable work 
of A. Priedmann [1,2], who solved the Einstein gravitational field equations and found that 
they admit non-static cosmological solutions presenting an expandeing universe. Einstein, 
believing that the universe should be static and unchanged forever, suggested a modification 
to his gravitational field equations by adding to them the so-called cosmological term which 
can stop the expansion. 

Soon after that E. Hubble [3,4] found experimentally that the far-away galaxies are 
receding from us, and that the farther the galaxy the bigger its velocity. In simple words, 
the universe is indeed expanding according to a simple physical law that gives the relationship 
between the receding velocity and the distance, 

V = HoR. (1) 

Equation (1) is usually referred to as the Hubble law, and Hq is called the Hubble constant. 
It is tacitly assumed that the velocity is proportional to the actual measurement of the 
redshift z of the receding objects by using the non-relativistic relation z = v/c, where c is 
the speed of light in vacuum. 

The Hubble law does not resemble standard dynamical physical laws that are familiar in 
physics. Rather, it is a cosmological equation of state of the kind one has in thermodynamics 
that relates the pressure, volume and temperature, pV — RT [5]. It is this Hubble's equation 
of state that will be extended so as to include gravity by use of the full Einstein theory of 
general relativity. The obtained results will be very simple, expressing distances in terms 
of redshifts; depending on the value of fl — p/pc we will have accelerating, constant and 
decelerating expansions, corresponding to Q < 1, Q = 1 and Q > 1, respectively. But the 
last two cases will be shown to be excluded on physical evidence, although the universe had 
decelerating and constant expansions before it reached its present accelerating expansion 
stage. As is well known, the standard theory does not deal with this problem. 

Before presenting our theory, and in order to fix the notation, we very briefly review the 
existing theory [6,7]. In the four-dimensional curved space-time describing the universe, our 
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spatial three-dimensional space is assumed to be isotropic and homogeneous. Co-moving 
coordinates, in which qqq — 1 and gok — 0, are employed [8,9]. Here, and throughout our 
paper, low-case Latin indices take the values 1, 2, 3 whereas Greek indices will take the 

values 0, 1, 2, 3, and the signature will be (-1 ). The four-dimensional space-time is 

split into 103 parts, and the line-element is subsequently written as ds^ — dt^ — dl^, where 
dl^ = ^^"^ Qkidx^dx^ = —gkidx^dx\ and the 3x3 tensor '^^^gu = —gu describes the geometry 
of the three-dimensional space at a given instant of time. In the above equations the speed 
of light c was taken as unity. 

Because of the isotropy and homogeneity of the three-geometry, it follows that the cur- 
vature tensor must have the form 

where X is a constant, the curvature of the three-dimensional space, which is related to the 
Ricci scalar by ^^^i? = —6K [10]. By simple geometrical arguments one then finds that 

dl^ = (l - r7i?2) dr^ + [dO'' + sin^ ^#2) , (3) 

where r < R. Furthermore, the curvature tensor corresponding to the metric (3) satisfies 
Eq. (2) with K = 1/B?. In the "spherical" coordinates {t,r,9,4>) we thus have 

gn = -{l-ryR'y\ (4) 

R is called the radius of curvature of the universe (or the expansion parameter) and its value 
is determined by the Einstein field equations. 

One then has three cases: (1) a universe with positive curvature for which K — l/i?^; 
(2) a universe with negative curvature, K = — 1/i?^; and (3) a universe with zero curvature, 
K — 0. The gii component for the negative-curvature universe is given by 

gu--{l + rVR'y\ (5) 

For the zero-curvature universe one lets — > 00. 
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Although general relativity theory asserts that all coordinate systems are equally vahd, 
in this theory one has to change variables in order to get the "right" solutions of the Einstein 
field equations according to the type of the universe. Accordingly, one makes the substitution 
r — Rsinx for the positive-curvature universe, and r — R sinh x for the negative-curvature 
universe. Not only that, but the time-like coordinate is also changed into another one rj by 
the transformation dt = Rdrj. The corresponding line elements then become: 

ds^ = R^ {t]) [d?]'^ - dx^ - sin^ x {dO^ + sin^ Odcj)^)] (6a) 

for the positive-curvature universe, 

ds^ = R^ {t]) [d7]^ - dx^ - sinh^ X (dO^ + sin^ ^^0^)] (66) 

for the negative-curvature universe, and 

ds^ = R^ {t]) [d?]"^ - dr^ - [dO'^ + sin^ Odcj)'^)] (6c) 

for the flat three-dimensional universe. In the sequel, we will see that the time-like coordinate 
in our theory will take one more different form. 

The Einstein field equations are then employed in order to determine the expansion 
parameter R{r)). In fact only one field equation is needed, 

ro _ 1^0^ + = SttGT", (7) 

where A is the cosmological constant, and c was taken as unity. In the Friedmann models 
one takes A = 0, and in the comoving coordinates used one easily finds that Tq — p, the 
mass density. While this choice of the energy-momentum tensor is acceptable in standard 
general relativity and in Newtonian gravity, we will argue in the sequel that it is not so for 
cosmology. At any rate, using p{t) — M/27r^i?^, where M is the "mass" and 27r^i?^ is the 
"volume" of the universe, one obtains 

3 [{dR/dtf + l] /R^ = AGM/nR^ + A, (8a) 
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or, in terms of 77 along with taking A = 0, 

3 [{dR/dr]f + B?] /R = AGM/t:, (9a) 
The solution of this equation is 

R^*R{l-cosri) , (10a) 

where *R is a constant, *R — 2G'M/37r, and from dt — Rdrj we obtain 

t ^*R{r) -sinrj) . (11a) 

Equations (10a) and (11a) are those of a cycloid, and give a full representation for the radius 
of curvature of the universe. 

Similarly, one obtains for the negative-curvature universe the analog to Eqs. (8a) and 
(9a), 

3 [{dR/dtf - 1] /R^ = AGM/ttR^ + A, (86) 
3 [{dR/dr)f - R"^] /R = AGM/n, (96) 
the solution of which is given by 

R = *R{coshr) -T)) , {10b) 

t = *R{smhr] -Tj) . (116) 

Finally, for the universe with a fiat three-dimensional space the Einstein field equations 
yield the analog to Eqs. (8a) and (9a), 

3 {dR/dtf lie = AGM/tvR^ + A, (8c) 

3 {dR/d7]f /R = AGM/n. (9c) 
As a function of t, the solution is 

R^{3GM/7rf^W/\ (10c) 



An extension of the Priedmann models was carried out by Lemaitre, who considered 
universes with zero energy- momentum but with a non-zero cosmological constant. One 
then again obtains positive-curvature, negative- curvature and zero- curvature (also called de 
Sitter) universes. While these models are of interest mathematically they have httle, if any, 
relation to the physical universe. 

In the final analysis, it follows that the expansion of the universe is determined by the 
so-called cosmological parameters. These are the mass density p, the Hubble constant H 
and the deceleration parameter q. We will not go through this, however, since we will con- 
centrate on the theory with dynamical variables that are actually measured by astronomers: 
distances, redshifts and the mass density. 

One of the Priedmann theory assumptions is that the type of the universe is determined 
by O = p/ Pc-i where Pc = 3Hq/8t:G, which requires that the sign of {fl — 1) must not change 
throughout the evolution of the universe so as to change the kind of the universe from one to 
another. That means in this theory, the universe has only one kind of curvature throughout 
its evolution and does not permit going from one curvature into another. In other words 
the universe has been and will be in only one form of expansion. It is not obvious, however, 
that this is indeed a valid assumption whether theoretically or experimentally. As will be 
shown in the sequel, the universe has actually three phases of expansion and it does go from 
one to the second and then to the third phase. 

A new outlook at the universe expansion can be achieved and is presented here. The new 
theory has the following features: (1) the dynamical variables arc Hubble's, i.e. distances and 
redshifts, the actually-measured quantities by astronomers; (2) it is fully general relativistic; 
(3) it includes two universal constants, the speed of light in vacuum c, and the Hubble 
time in the absence of gravity r (might also be called the Hubble time in vacuum); (4) the 
redshift parameter z is taken as the time-like coordinate; (5) the energy-momentum tensor 
is represented differently; and (6) it predicts that the universe has three phases of expansion: 
decelerating, constant and accelerating, but it is now in the stage of accelerating expansion 
phase after passing the other two phases. 
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Our starting point is Bubble's cosmological equation of state, Eq. (1). One can keep the 
velocity v in equation (1) or replace it with the redshift parameter z by means oi z — v/c. 
Since R= {xi,X2,X3), the square of Eq. (1) then yields 

c'H^h''-(xl + xl + xl)=Q. (12) 

Our aim is to write our equations in an invariant way so as to enable us to extend them to 
curved space. Equation (12) is not invariant since H^^ is the Hubble time at present. At 
the limit of zero gravity, Eq. (12) will have the form 

cV^z^ - (xj + xl + xf) = 0, (13) 

where r is Bubble's time in vacuum, which is a universal constant the numerical value 
of which will be determined in the sequel by relating it to Hq^ at different situations. 
Equation (13) provides the basis of a cosmological special relativity and has been investigated 
extensively [11-16]. 

In order to make Eq. (13) adaptable to curved space we write it in a differential form: 

c^r'^dz^ - [dxl + dxl + dxfj = 0, (14) 
or, in covariant form in flat space, 

ds^ = rj^^dx^dx^ = 0, (15a) 

where rj^i, is the ordinary Minkowskian metric, and our coordinates arc {x^ , x^ , x'^ , x^) = 
{crz,Xi,X2,X3). Equation (15a) expresses the null condition, familiar from light propaga- 
tion in space, but here it expresses the universe expansion in space. The generalization of 
Eq. (15a) to a covariant form in curved space can immediately be made by replacing the 
Minkowskian metric rj^i, by the curved Riemannian geometrical metric g^,^, 

ds"^ = g^.^dx^'dx" = 0, (156) 

obtained from solving the Einstein field equations. 

7 



Because of the spherical symmetry nature of the universe, the metric we seek is of the 
form [8] 

ds^ = c^^dz" - e^dr^ - (dO^ + sin^ Odcp") , (16) 

where co-moving coordinates, as in the Priedmann theory, are used and A is a function of 
the radial distance r. The metric (16) is static and solves the Einstein field equation (7). 
When looking for static solutions, Eq. (7) can also be written as 

e-^ (A'/r - 1 /r') + l/r^ = SttGT^ ( 17) 

when A is taken zero, and where a prime denotes differentiation with respect to r. 

In general relativity theory one takes for Tq = p. So is the situation in Newtonian gravity 
where one has the Poisson equation = inGp. At points where p = one solves the 
vacuum Einstein field equations and the Laplace equation = in Newtonian gravity. 
In both theories a null (zero) solution is allowed as a trivial case. In cosmology, however, 
there exists no situation at which p can be zero because the universe is filled with matter. 
In order to be able to have zero on the right-hand side of Eq. (17) we take not as equal 
to p but to p — Pci where Pc is chosen by us now as a constant given by Pc — S/SttGt^. This 
approach has been presented and used in earlier work [17]. 

The solution of Eq. (17), with Tq — p — pc, is given by 

e-^ = 1 - (Q - 1) rVcV^, (18) 

where ^ — pjpc- Accordingly, if Q > 1 we have grr — — (1 — ^^/-R^)~^, where 

i?2 ^ (^^ - 1) , (19a) 

exactly equals to gw given by Eq. (4) for the positive-curvature Priedmann universe that 
is obtained in the standard theory by purely geometrical manipulations. If Q < 1, we can 
write gf^r = — (1 + r"^ /B?)~^ with 

= c^T^I (1 - 1^) , (196) 
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which is equal to gu given by Eq. (5) for the negative-curvature Priedmann universe. In 
the above equations r < R. 

Moreover, we know that the Einstein field equations for these cases are given by Eqs. 
(8) which, in our new notation, have the form 



As is seen from these equations, if one neglects the first term in the square brackets with 
respect to the second ones, will be exactly reduced to their values given by Eqs. (19). 

The expansion of the universe can now be determined from the null condition ds — 0, 
Eq. (15b), using the metric (16). Since the expansion is radial, one has dO = dcf) = 0, and 
the equation obtained is 



The second term in the square brackets of Eq. (21) represents the deviation from constant 
expansion due to gravity. For without this term, Eq. (21) reduces to dr/dz — ct or 
dr/dv = T, thus r = rv+constant. The constant can be taken zero if one assumes, as usual, 
that at r = the velocity should also vanish. Accordingly we have r = tv or v = T~^r. 
When D, — 1, namely when p — pc, we have a constant expansion. 

The equation of motion (21) can be integrated exactly by the substitutions 




(20a) 




(206) 



dr/dz ^CT l + (l-n)rVcV2 



n 1/2 



(21) 



sinx = ar/cT; fl > 1, 



(22a) 



sinhx = Pr/cr; fl <1, 



(226) 



where 



a — 



(Q-l)^/^ /3=(i-o) 



1/2 



(23) 



For the Q > 1 case a straightforward calculation, using Eq. (22a), gives 



dr — {cr / a) cos xdx, 



(24) 
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and the equation of the universe expansion (21) yields 

dx — adz. (25a) 

The integration of this equation gives 

X — CKZ + constant. (26a) 

The constant can be determined, using Eq. (22a). For at x = 0, we have r — and z — 0, 
thus 

X = az, (27a) 
or, in terms of the distance, using (22a) again, 

r (z) = (cr/a) sin az; a ^ (Q - 1)^^^ . (28a) 

This is obviously a decelerating expansion. 

For Q < 1, using Eq. (22b), then a similar calculation yields for the universe expansion 
(21) 

dx = Pdz, (256) 

thus 

X — Pz + constant. (266) 
Using the same initial conditions used above then give 

X = f3z, {27b) 

and in terms of distances, 

r{z) = {cT/p)smh(3z; /3=(l-0)^/^ (286) 

This is now an accelerating expansion. 
For n = 1 we have, from Eq. (21), 

d'^r/dz'^ = 0. (25c) 
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The solution is, of course, 



r {z) 



= CTZ. 



(28c) 



This is a constant expansion. 

It will be noted that the last solution can also be obtained directly from the previous two 
ones for Q > 1 and ^2 < 1 by going to the limit 2; — > 0, using L'Hospital lemma, showing that 
our solutions are consistent. It will be shown later on that the constant expansion is just 
a transition stage between the decelerating and the accelerating expansions as the universe 
evolves toward its present situation. 

Figure 1 describes the Hubble diagram of the above solutions for the three types of 
expansion for values of from 100 to 0.24. The figure describes the three-phase evolution 
of the universe. Curves (1) to (5) represent the stages of decelerating expansion according 
to Eq. (28a). As the density of matter p decreases, so does fl also, the universe goes over 
from the lower curves to the upper ones, and it does not have enough time to close up to a 
big crunch. The universe, subsequently, goes to curve (6) with Q = 1, at which time it has 
a constant expansion for a fraction of a second. This then followed by going to the upper 
curves (7) and (8) with Q < 1 where the universe expands with acceleration according to 
Eq. (28b). A curve of this kind fits the present situation of the universe. 

In order to decide which of the three cases is the appropriate one at the present time, we 
have to write the solutions (28) in the ordinary Hubble law form v — H^r. To this end we 
change variables from the redshift parameter z to the velocity v by means oi z = v/c for v 
much smaller than c. For higher velocities this relation is not accurate and one has to use a 
Lorentz transformation in order to relate z io v. A simple calculation then shows that, for 
receding objects, one has the relations 



z — 



[{l+v/c)/{l-v/ct''-l, 



(29a) 




(296) 



We will assume that v <^c and consequently Eqs. (28) have the forms 



r {v) — {ct/o) sin {av/c) , 



(30a) 
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r (v) — {ct/ (5) sinh {(5v/c) , 
r {v) = TV. 



(306) 
(30c) 



Expanding now Eqs. (30a) and (30b) and keeping the appropriate terms, then yields 



for the Q > 1 case, and 



r — TV (l — a^v^ , 



TV (l + /3V/6c2), 



(31a) 



(316) 



for Q < 1. Using now the expressions for a and (3, given by Eq. (23), in Eqs. (31) then 
both of the latter can be reduced into a single equation 



r — TV 



1 + (1 - O) v'^/Qc^ 



(32) 



Inverting now this equation by writing it in the form v = Hqt, we obtain in the lowest 
approximation for Hq the following: 



HQ^h 1 - (1 - Q) ^;V6c^ 



where h — t . Using v t/t, or z — v/c, we also obtain 



(33) 



Ho^h l-(l-Q)rV6cV 



= h 



1 - (1 - Q) zV6 



(34) 



Consequently, Hq depends on the distance, or equivalently, on the redshift. As is seen, 
Hq has meaning only for r — >^ or 2; — > 0, namely when measured locally, in which case it 
becomes h. 

In recent years observers have argued for values of Hq as low as 50 and as high as 90 
km/sec-Mpc, some of the recent ones show 80 ± 17 km/sec- Mpc [18-26]. There are the 
so-called "short" and "long" distance scales, with the higher and the lower values for Hq 
respectively [27]. Indications are that the longer the distance of measurement the smaller 
the value of Hq. By Eqs. (33) and (34) this is possible only for the case in which < 1, 
namely when the universe is at an accelerating expansion. 
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Figures 2 and 3 show the Hubble diagrams of the predicted by theory distance-redshift 
relationship for the accelerating expanding universe at present time, whereas figures 4 and 
5 show the experimental results [28-29]. 

Our estimate for h, based on pubhshed data, is /i 85 — 90 km/sec- Mpc. Assuming 
^ 85 km/sec-Mpc, Eq. (34) then gives 

Ho^h[l-1.3xl0-^(l-n)r^], (35) 

where r is in Mpc. A computer best-fit can then fix both h and Q. 

In summary, we have presented a general-relativistic theory of cosmology, the dynamical 
variables of which are those of Bubble's, i.e. distances and redshifts. The theory describes 
the universe as having a three-phase evolution with a decelerating expansion, followed by 
a constant and an accelerating expansion, and it predicts that the universe is now in the 
latter phase. As the density of matter decreases, while the universe is at the decelerating 
phase, it does not have enough time to close up to a big crunch. Rather, it goes to the 
constant expansion phase, and then to the accelerating stage. The equations obtained for 
the universe expansion are very simple. 

Finally, it is worth mentioning that if one assumes that the matter density at the present 
time is such that Qq = 0.24, then it can be shown that the time at which the universe 
had gone over from a decelerating to an accelerating expansion, i.e. the constant expansion 
phase, occured at cosmic time 0.03t. 

The idea to express cosmological theory in terms of directly-measurable quantities, such 
as distances and redshifts, was partially inspired by Albert Einstein's favourite remarks on 
the theory of thermodynamics in his Autobiographical Notes [30] . 
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FIGURE CAPTIONS 



Fig. 1 Hubble's diagram describing the three-phase evolution of the universe according to 
Einstein's general relativity theory. Curves (1) to (5) represent the stages of decelerating 
expansion according to r{z) = (cr/a) sinaz, where a = {Q. — 1)"^^^, ^ = p/Pc: with pc a 
constant, pc — S/SttGt^, and c and r are the speed of light and the Hubble time in vacuum 
(both universal constants). As the density of matter p decreases, the universe goes over 
from the lower curves to the upper ones, but it does not have enough time to close up to 
a big crunch. The universe subsequently goes to curve (6) with D, — 1, at which time it 
has a constant expansion for a fraction of a second. This then followed by going to the 
upper curves (7)-(8) with fl < 1 where the universe expands with acceleration according 
to r{z) — (ct/P) sinhpz, where P — {1 — One of these last curves fits the present 

situation of the universe. 

Fig. 2 Hubble's diagram of the universe at the present phase of evolution with accelerating 
expansion. 

Fig. 3 Hubble's diagram describing decelerating, constant and accelerating expansions in a 
logarithmic scale. 

Fig. 4 Distance vs. redshift diagram showing the deviation from a constant toward an 
accelerating expansion. [Sourse: A. Riess et al., Astron. J. 116, 1009 (1998)]. 

Fig. 5 Relative intensity of light and relative distance vs. redshift. [Sourse: A. Riess et al., 
Astron. J. 116, 1009 (1998)]. 
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Fig. 4. 
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